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Abstract 

The transition process induced by the 
interaction of an isolated roughness with 
acoustic disturbances in the free stream is 
numerically investigated for a boundary layer 
over a flat plate with a blunted leading edge at a 
free stream Mach number of 3.5. The roughness 
is assumed to be of Gaussian shape and the 
acoustic disturbances are introduced as boundary 
condition at the outer field. The governing 
equations are solved using the 5 th -order accurate 
weighted essentially non-oscillatory (WENO) 
scheme for space discretization and using third- 
order total- variation-diminishing (TVD) Runge- 
Kutta scheme for time integration. The steady 
field induced by the two and three-dimensional 
roughness is also computed. The flow field 
induced by two-dimensional roughness exhibits 
different characteristics depending on the 
roughness heights. At small roughness heights 
the flow passes smoothly over the roughness, at 
moderate heights the flow separates downstream 
of the roughness and at larger roughness heights 
the flow separates upstream and downstream of 
the roughness. Computations also show that 
disturbances inside the boundary layer is due to 
the direct interaction of the acoustic waves and 
isolated roughness plays a minor role in 
generating instability waves. 

Introduction 

Transition from laminar to turbulent state in 
shear flows occurs due to evolution and 
interaction of different disturbances inside the 
shear layer. Though there are several 
mechanisms and routes to go from a laminar to a 
turbulent state, all of them generally follow 
these fundamental processes: 

• Receptivity 


• Linear instability 

• Nonlinear instability and saturation 

• Secondary instability and breakdown 
to turbulence 


In the receptivity process, the unsteady 
disturbances in the environments such as 
acoustic and turbulence interact with the 
inhomogeneities in the geometry such as 
roughness and generate instability waves inside 
the shear layer. In quiet environments, the initial 
amplitudes of these instability waves are small 
compared to any characteristic velocity and 
length scales in the flow. In the second stage, the 
amplitudes of these instability waves grow 
exponentially downstream and this process is 
governed by the linearized Navier-Stokes 
equation. Further downstream, the amplitudes of 
the disturbances become large and the nonlinear 
effects inhibit the exponential growth and the 
amplitude of the waves eventually saturate. In 
the next stage, this finite amplitude saturated 
disturbances become unstable to two- and/or 
three-dimensional disturbances. This is called 
secondary instability and beyond this stage the 
spectrum broadens, due to complex interactions 
and further instabilities, and the flow becomes 
turbulent in a short distance downstream. In this 
paper, the receptivity process induced by an 
isolated roughness and acoustic disturbances in a 
supersonic boundary layer with a free stream 
Mach number of 3.5 is investigated. 

Predicting transition onset and transition end 
points accurately, modeling this transitional 
region and modeling the turbulence region are 
the major difficulties in computing the 
aerodynamic quantities accurately using the 
CFD codes. Our understanding of different 
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instability mechanisms and of different 
transition processes in shear layers have greatly 
improved in the last several decades. However 
the transition prediction methods have not made 
much progress. The main difficulty is due to the 
nature of the transition process itself. The 
transition process mainly depends on the 
boundary layer characteristics and on the 
frequency and wave number distributions of the 
disturbances that enter the boundary layer. The 
laminar boundary layer profiles can be 
computed easily. The problem is computing, 
predicting or prescribing the initial spectral, 
amplitude and phase, distribution of the 
disturbances inside the boundary layer. As we 
discussed earlier, the initial disturbances are 
generated by the interaction of the free stream 
unsteady disturbances and the roughness on the 
surface. These two are stochastic in nature and 
are difficult to quantify in general. In any new 
transition prediction strategy, one should 
quantify these two quantities and should 
determine what is the minimum amount of 
information necessary to predict the transition 
onset accurately. The objectives of this research 
work are to answer some of these questions and 
eventually to come up with an improved 
transition prediction method. 

There are numerous investigations conducted on 
the interaction of acoustic waves with 
supersonic boundary layers. Mack (1976), 
Gapanov (1977) investigated the interaction of 
acoustic waves with a supersonic boundary layer 
at finite incident angles using inhomogeneous 
stability equations. One important finding was 
that due to the interaction, the acoustic waves 
excite disturbances inside the boundary layer, 
which are much larger than that in the free 
stream. Gapanov and Smorodsky (1999) studied 
the interaction of streamwise acoustic waves 
with a non-parallel boundary layer. The analysis 
and the calculations showed that the 
disturbances inside the boundary layer reach 
significant values compared to that in the 
outside. It was also observed that there exists a 
critical Reynolds number where this excitation is 
the highest. Fedorov and Khokhlov (1991) 
investigated using asymptotic theory the 
excitation of first and second modes by the 
acoustic waves near the leading edge region. 


Gapanov (1993) showed that the energy is 
transmitted to the T-S waves near the critical 
region where the ratio between the reflected 
waves and the incident wave is the largest. There 
are not many numerical calculations performed 
to investigate the interaction of acoustic waves 
with a flat plate including bluntness. Since there 
are forced disturbances and the modal 
disturbances what are the relative magnitudes of 
them and how far the forced disturbances persist 
are not very well understood. 

There were several transition experiments 
performed at NASA Langley in the Mach 3.5 
Supersonic Low-Disturbance tunnel. Boundary- 
layer transition data on a flat plate and on a 
cone, and free stream noise levels and the power 
spectral distribution of the free stream noise are 
presented in Chen et. al (1988). As a first step, 
the computations are performed for the same 
conditions as in the experiment. An isolated two 
and three-dimensional Gaussian shaped 
roughness is placed near the neutral stability 
region and the steady flow field induced by them 
are investigated. Secondly the interaction of the 
two dimensional acoustic disturbance with a 
blunted semi-infinite plate with and without 
two-dimensional roughness element located near 
the neutral point are investigated to identify the 
effect of acoustic incident angle and the 
roughness in receptivity process in a supersonic 
boundary layer. A schematic diagram of the 
computational set up is depicted in figure 1. 


Governing Equations. 


The equations solved are the three-dimensional 
unsteady compressible Navier-Stokes equations 
in conservation form 
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Here (x,y,z) are the Cartesian coordinates, (u ,v 
,w) are the velocity components, p is the density, 
and p is the pressure. E is the total energy given 
by 


E- e + ■ 


u 2 + v 2 + w 2 


e = c v T , p=pRT. 


Here e is the internal energy and T is the 
temperature. The shear stress and the heat flux 
are given by 



The viscosity (p) is computed using 
Sutherland’s law and the coefficient of 
conductivity (k) is given in terms of the Prandtl 
number Pr. The variables p, p, T and velocity 
are non-dimensionalised by their correspo nding 
reference variables p«, p a , T« and JRT X 
respectively. The reference value for length is 
computed by -Jvx 0 TeT , where x 0 is the 

location of the beginning of the computational 
domain from the leading edge in the streamwise 
direction. For the computation, the equations are 
transformed from physical coordinate system (x, 
y, z) to the computational curvilinear coordinate 
system (g ,r ]£ ) in a conservative manner and 
the governing equations become 



The components of the flux in the computational 
domain are related to the flux in the Cartesian 
domain by 



where J = 


r <?(g,q,g) 1 
[<?(jr,y,z) J' 


Solution Algorithm 

The governing equations are solved using the 
5' h -order accurate weighted essentially non- 
oscillatory (WENO) scheme for space 
discretization and using third-order total- 
variation-diminishing (TVD) Runge-Kutta 
scheme for time integration. The WENO and the 
TVD methods and the formulas are explained in 
Shu (1992) and the application of ENO method 
to N-S equations is given in Atkins (1991). The 
solution method used is described in an earlier 
paper Balakumar et. al. (2002). 

The height distribution of a rough surface is 
generally determined using a Gaussian 
probability distribution and a correlation length 
(Saillard and Sentenac 2001). Before a general 
rough surface is considered, in this paper the 
disturbances induced by isolated roughness are 
studied. The roughness is assumed to be of 
Gaussian shape in the form 

y - f(x,z) 


where h is the height of the roughness, cr, and 
o 2 determine the width of the roughness in the 
streamwise and spanwise directions and x=x 0 is 
the location of the roughness. The grid is 
generated by transforming the surface into a 
straight line using the expression 

y-f(x,z) 

H(x) - f(x,z) ’ 

where y = H(x) is the outer boundary in the 
computational domain. 
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The acoustic field that impinges on the outer 
boundary is taken to be in the following form. 

p = Re al{ pe^iPz^y-iax j_ 

Here a ac , (> ac . e ac are the acoustic wavenumber, 
and co is the frequency of the acoustic 
disturbance. 

Results 

First the flow field induced by two and three- 
dimensional isolated Gaussian shaped roughness 
elements of different heights are computed. 
Secondly, the disturbances generated by the 
interaction of two-dimensional acoustic 
disturbances with a blunted flat plate with and 
without roughness are computed. The flow 
parameters are given in Table 1. 


Non-dimensional frequency F = 1*10' 3 =41.0 

kHz 

The non-dimensional frequency F is defined as 

vl 

where f is the frequency in Hertz. 

The grid is generated using analytical conformal 
mapping formulae. The grid is clustered in the r| 
directions close to the wall and near the critical 
layer region and is also clustered in the § and £ 
directions close to the roughness and is shown in 
figure 2. We present the results for the 
parameters x 0 =1 in., a, = o 2 = .25. here x 0 is 
the location of the roughness element from the 
leading edge. 

Linear instability 


Table 1 . 

Flow parameters for the wind tunnel model . 

Freestream Mach number M, = 3.5 

Freestream Reynolds number Re a = 12*10 6 /ft 

Freestream density p„ = 2.249* 10' 2 lbm/ft 3 

Freestream pressure p, = 187.74 lbf/ft 2 

Freestream velocity U« = 2145.89 ft/s 

Freestream temperature T x = 156.42 °R 

Free stream kinematic viscosity v x = 1.7882*1 O' 4 * 
ft 2 / s 

Wall temperature = Adiabatic condition 
Prandtl number Pr = 0.72 
Ratio of specific heats y = 1 .4 

Length scale ^^ = 5.892*10 ^ ft. 
(x 0 = 0.5 in.) 

The boundary layer thickness at X=1 in. 8 0 = 
.01275 in. 


As a prelude for future reference, in figure 3 the 
linear stability results for the boundary layer 
over a flat plate is presented. The figure depicts 
the neutral stability diagram in (Re, F), (Re, a) 
and (Re, (3) planes for different wave angles 0, 
45, 60 and 70 degrees. The figure also shows the 
N-Factor curves and the growth rates for the 
most amplified disturbances. The critical 
Reynolds number is about 193 and this occurs 
for an oblique wave of angle 60 degrees. The 
most amplified frequency is about F= 1.0-1. 25e- 
5 and the most amplified wave has a spanwise 
wavenumber of p=. 025 and this corresponds to 
about .178 inches in dimensional units and is 
equivalent to about 14 boundary layer 
thicknesses. It is also observed that at higher 
Reynolds numbers Re > 1000, only the low 
frequency disturbances F < 3.E-5 are unstable. 
This implies that acoustic disturbance with 
frequencies less than 120 kHz may be the 
relevant frequency range for generating 
instability waves inside the boundary layer. The 
frequency of the most amplified wave is about 
40 kHz. The maximum N-factor at X=12 in, 
(Re = 3464) is about 8.6. 

Steady flow due to the roughness. 
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In these computations, an isolated two and three- 
dimensional roughness are placed at Xo=l in. 
from the leading edge. The Reynolds number at 

this location is ^' xA ° - = 1000. and the 


v„ 


parameters in the Gaussian 
distribution a, = 0 2 = 0.25 . This gives the 
width of the roughness element as half the 
boundary layer thickness. 


Two-dimensional roughness. 


Figures 4-8 show the results for the two- 
dimensional roughness elements. The 
computations are performed for different 
roughness heights h/8 0 =l/16, 1/8, 1/4 and 1/2. 
Figure 4 shows the contours of density inside the 
boundary layer and in the inviscid region for two 
cases h/6 = 1/16 and 1/2. Steady solutions are 
obtained in all the cases and there is no vortex 
shedding observed. For h/6 = 1/2, in the inviscid 
region the results clearly show a compression 
wave in the front, an expansion wave in the 
middle and another compression wave at the 
end. As expected these waves are weaker for the 
small roughness h/8=l/16 compared to h/8=l/2. 
Figure 5 shows the contours of the normal 
velocity (V) and the streamline patterns inside 
the boundary layer near the roughness for 
roughness heights h/6 = 1/16, 1/8, 1/4 and 1/2. 
The interesting observation is the type of 
separation induced by the roughness. At small 
roughness heights h/8 = 1/16 there is no 
separation observed and the flow goes over the 
roughness smoothly. For roughness heights 1/16 
< h/8 < 1/4, flow separates downstream of the 
roughness. The flow separates slightly 
downstream of the peak of the roughness and for 
h/8 = 1/8 it reattaches at about 4 roughness 
heights downstream of the peak. For roughness 
heights h/8 > 1/4, the flow separates in the 
upstream and downstream of the roughness. The 
length of the separated regions increases with 
roughness height and for h/8 = 1/2, the upstream 
separated region is about 1 5 h and the 
downstream region is about 8h. 


Figure 6 shows the variation of the normal 
density in the streamwise direction (x) at 
constant heights y/8 0 = 1 and 2 for different 
roughness heights h/8 0 = 1/16, 1/8, 1/4, 1/2. As 
observed earlier, a smooth compression in the 
upstream, a sharp and a strong expansion fan 
across the roughness peak and a strong 
compression downstream are observed. The 
strength of these waves increases with 
increasing heights. In figures 7 and 8, the 
velocity profiles at different streamwise 
locations relative to the roughness location are 
presented and they are compared with the 
unperturbed Blasius velocity profiles. Figure 7 
shows the difference between the velocity 
profiles induced by the roughness and the 
Blasius profiles at the upstream streamwise 
locations (x-x 0 )/8 0 = -16, -4, -2 and at the 
downstream locations (x-x 0 )/8o = 2, 4, 16. The 
first observation is that the modification to the 
boundary layer is more in the region upstream of 
the roughness than in the downstream of the 
region. The change in maximum velocity is 
about 1 % at the upstream location (x-x 0 )/8„ = -2 
and is about 0.3% at the downstream location (x- 
x 0 )/8 0 = 2. The shape of the modification to the 
velocity profiles is also different from the 
upstream to the downstream. In the downstream, 
the velocity is reduced more near the wall and 
near the outer part of the boundary layer than in 
the middle of the boundary layer and this may 
have a strong influence in the stability 
characteristics. Figure 8 shows the velocity 
profiles and the comparison with the Blasius 
profiles for the roughness h/8 0 = 1/2. Again the 
modification is larger in the upstream region of 
the roughness about 40% at (x-x 0 )/8o = -4 
compared to that in the downstream region 
which is about 15% at (x-x 0 )/8 0 = 4. The 
difference from the Blasius profiles decreases 
with increasing distance from the roughness and 
it is less than 1 % at (x-x 0 )/8 0 = 20. 

Three-dimensional roughness. 

Figures 9 and 10 show the flow field induced by 
three-dimensional roughness of heights h/8 = 
1/16 and 1/8. In the spanwise direction periodic 
boundary condition is applied. The wavenumber 
in the spanwise direction is p=.025, which 
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correspond to the most, amplified disturbance. 
Figure 9 shows the streamline patterns and the 
contours of the density in the (x-y) plane at z=0. 
In these cases also, steady solutions are obtained 
and no vortex shedding is observed. For the 
roughness height h/6 0 =l/8 a small recirculation 
region is observed downstream. In the three- 
dimensional case the compression and the 
expansion waves are weaker compared to the 
two-dimensional case. Figure 10 shows the 
contours of the spanwise velocity component 
(w) in the (y-z) planes at different streamwise 
locations relative to the roughness location, (x- 
x 0 )/S 0 =-2, 0., 7.5 and 190. for the roughness 
height h/6 0 =l/16. Figures 10(a-c) show the flow 
field near the roughness height and the figure 
lOd shows the results far downstream from the 
roughness. The maximum spanwise velocity is 
about .02 and within about two boundary layer 
thicknesses upstream and ten boundary layer 
thicknesses downstream the velocities are 
reduced by two orders of magnitude. Figure lOd 
shows the flow field induced by the roughness in 
the inviscid part of the flow. Due to the periodic 
and symmetric boundary conditions at the 
boundaries the Mach waves radiated in the 
spanwise directions bounce back and forth 
inside the domain and this is clearly seen in 
figure lOd in the outer part of the flow field. 
This made the convergent to the steady state 
very slow. 


Acoustic waves. 


The linearized Euler equations in a uniform 
mean flow are: 


dp rT dp du dv dw ^ 
— + t/ 0 — + p° — + ft — + p° — = 0, 

dt dx dx dy dz 

du , _ du dp 

Po + AM) T” = “ 

dt dx dx 

dv tj dv dp 

Pot +PoUo— =~* 
dt dx cry 

dw TJ dw dp 

Po — + Pof / o7 = -— - 
dt dx dz 


cJT 


dr dp 


A c- + p 0 f/ 0 c -=- +( ; 0 


dp 
dx ’ 


P o = Po*T 0 , 
p = p 0 RT + pRT 0 . 

The solution of this system can be written as 


a o 

Of 


PoW-") 

Pac 


Po(a ac U 0 -w) 

£ 


PoiVacUo-M) 

(r-i)7 0 


Po a o 

Here the pressure p is in the form 


P- 


P = Pamp e 


i{a ac x + fi M y+e ac z-u>t) 


The dispersion relation among the wavenumbers 
a ac , P ac , e ac an d the frequency oo is given by 


(cx ac U 0 - (of = (aj + + ej)al . 


For zero sweep acoustic disturbances £ a c=Ch the 
wavenumber a ac can be expressed as 


co cos 0 


or = - 


(U 0 cosd v ±a 0 ) 
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Here 0 v 


to \ 


tan' 


A 


\ a oc/ 


is the incident angle and 


for p ac < 0 the plate is radiated from above and 
for p ac > 0 the plate is radiated from below or 
represents the wave which is reflected from the 
plate (figure 11). The plus sign corresponds to 
the fast moving wave and the minus sign 
corresponds to the slow moving wave. The 
corresponding phase speeds are 

C = U n ± — — — . The wavenumber of the fast 
0 cos 0 V 


CO 

moving wave is a nr < 

U 0 + a 0 

(jo 

moving wave a ac > 

t/o ~ a o 

angle is limited by 6 y < cos 1 


- and for the slow 

and the incident 

(—). At M=3.5, 
\MJ 


F=1.25E-5 the wavenumbers of the fast moving 
and slow moving waves are a ac < 0.00687 and 
a ac > 0.01237 and the incidence angle of the 
slow moving wave is limited to 73.39 degrees. 
The wavenumber of the acoustic disturbances 
a ac for different incident angle 0 y is given in 
Table 2 for two different spanwise wavenumbers 


Table 2. 

Values of ou for different inclination ang le 


P=0.025 instability waves are plotted in figure 
12. The wavenumber for the two-dimensional 
wave decreases from .0124 at a Reynolds 
number of 220 to .01 19 at a Reynolds number of 
2000. For the three-dimensional disturbance it 
decreases from .01775 at a Reynolds number of 
190 to .01285 at a Reynolds number of 2000. 
The Reynolds number at the neutral points is 
820 and 610 respectively for the two and three- 
dimensional disturbances and they are also 
marked in figure 12. In this figure the range of 
the acoustic wavenumber a ac for the two and 
three-dimensional acoustic waves are also 
depicted. It is seen that perfect matching in the 
wavenumbers occur at very low Reynolds 
number range where the disturbances are 
marginally stable. For the unit Reynolds number 
of 12* 10 6 , Reynolds number of 200 is located at 
.04 inches from the leading edge of the plate. 
Hence there may be strong generation of 
instability waves in the nose region where the 
bluntness effects will be important. In the 
neutral and unstable region the acoustic and the 
instability wavenumbers do not match and some 
tuning must occur due to nonparallel effects of 
the boundary layer or due to some 
inhomogenities on the surface such as rough 
nesses. In this paper the interaction of the 
acoustic waves with the boundary layer with and 
without roughness in the neutral and unstable 
region is investigated. 


0 y 

E ac =0.0 

e ac =.025 

0 

.01237 

.01757 

5 

.01239 

.01758 

10 

.01245 

.01762 

20 

.01270 

.01780 

30 

.01319 

.01816 

40 

.01409 

.01886 


The variation of the wavenumber awith the 

U X 

Reynolds number |— — for a two- 


dimensional P=0 and a three-dimensional 


Acoustic waves and the continuous spectrum 

The unsteady flow field in a boundary layer 
consists of all the discrete eigensolutions and the 
continuous spectrum, which appears because the 
domain is unbounded and the linearized stability 
equations admit solutions that are bounded at 
infinity. The details about the derivation and the 
computation of the discrete and the continuous 
spectrum are given in Balakumar and Malik 
(1992). For a supersonic boundary layer, there 
exist seven branches of the continuous spectrum 
in the complex wavenumber space. Two of them 
are the fast and slow moving acoustic waves. 
Hence the free-stream acoustic disturbances 
enter the boundary layer through the continuous 
spectrum. 
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The solution of the linearized stability equations 
in the free stream can be written as 

where q i (i=l,8) are column vectors and 


- fj. 2 = a 2 + fi 2 - (a - (o) 2 M 2 . 

This is same as the dispersion relation for the 
acoustic waves and the inviscid continuous 
spectrum for a two-dimensional disturbance P=0 
is given by 


r, = r 2 = -p 2 , r 3 = -p 3 , r 4 = -p 4 , 
r 5i =&, h=p 2 , r 7 =p 3 , r g =p 4 , 

ft] = A] , = A 2 , /^ 3 A 3 , p 4 A 4 , 

A, = A, = i(e - (jo ) Re+ a 2 + /3 2 , 



+ /> 


23 



+ b 33 ) + 4(/?23^32 



2 


£>22 


+ ^23 


3 

— [(^22 + ^ 33 ) + 4 (^ 23^32 — ^ 22 ^ 33)] 2 > 


to 


Mcos0 v 


The eigenfunction corresponding to the 
continuous spectrum is obtained by solving the 
linearized stability equations with an 
inhomogeneous boundary condition at the far 
field. Keeping only the bounded solution at the 
far field and discarding the exponentially 
growing solution, the solution in the far field 
corresponding to the acoustic wave continuous 
spectrum takes the form 


b 22 =a 2 + /3 2 -Re(a -to) 2 

—yM 2 -(y-l)Pr M 2 

4 

3 9 ’ 

— Re+ i(a-co)yM 2 

Re(a - to) 2 j-^ - P f | 

^23 = 3 ’ 

— Re+ i(a-co)yM 2 

4 

b 32 = -i(a ~(o)(y - l)Pr A/ 2 Re, 
b 33 = i(a - to)PrRe+a 2 + /3 2 . 

The continuous spectrum in the a -plane are 
determined by solving these equations with 

A, = — A| 0 , A 3 = -A 30 , A 4 = — A^, 

where 0 <. A i0 , A 30 , A^, <. 00 . The continuous 
spectrum correspond to acoustic waves are 
obtained from 

A 4 = — A40 — —f-i . 

In the limit of large Reynolds number this 
equation takes the form 


u 

v 


w 

T 

P 

u 

w' 

T I 


= Cm** + C 2 q 2 e^ + C 3 q ^ + C 4 q 4 e m + C,q/ m 


The first three terms in the right hand side 
represent the exponentially decaying solution 
and the last two terms represent the oscillatory 
acoustic disturbances. Between the last two 
terms, the first term is the incident acoustic wave 
and the second term is the reflected acoustic 
wave. Hence by fixing the amplitude of the 
incident wave, for instance the amplitude of the 
pressure of the incoming wave is selected as 
C 4 q 45 = Pq, the coefficient C 4 can be fixed. 
Since there are eight equations with four 
unknowns Q, C 2 , C 3 and C 5 , eliminating C,, C 2 , 
C 3 and C 5 the necessary four inhomogeneous 
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boundary conditions at the far field are obtained. 
After the eigenfunction is solved for, the values 
of C„ C 2 , C 3 and C 5 can be calculated, 
especially the amplitude and the phase of the 
reflected acoustic wave C 5 can be inferred. 

Figure 1 3 shows the amplitude of the coefficient 
of the reflected wave C 5 as a function of the 
incidence angle at different Reynolds numbers 
for a two-dimensional and three-dimensional 
acoustic wave. For a two-dimensional wave it is 
seen that at low Reynolds numbers the reflection 
coefficient is very large at small incident angles. 
It is about 27 at 3 degrees at a Reynolds number 
of 100 and it decreases with increasing Reynolds 
numbers. For three-dimensional acoustic 
disturbances the maximum reflection coefficient 
occurs at larger incident angles compared to the 
two dimensional case. In figure 14, the 
amplitude of the density fluctuations obtained by 
solving the inhomogeneous problem with the 
acoustic forcing at the free stream is presented. 
The results show the distribution for two 
incident angles 1.5 and 10 degrees at two axial 
locations 0.1 and 4.0 inches from the leading 
edge. The first observation is that with 
increasing Reynolds number the direct acoustic 
response is decreasing and the second is that at 
the larger incident angle of 10 deg. the 
disturbance inside the boundary layer is an order 
of magnitude larger than at the small incident 
angle of 1.5 deg. Since the low Reynolds 
numbers occur very close to the leading edge, 
the bluntness effect and the non-parallel effect 
will influence these results. 

As it is discussed in the introduction, the 
instability waves are generated inside the 
boundary layer by the free-stream disturbances 
and the roughness on the surface. Since the 
acoustic disturbances have the wavelengths in 
the same range as the unstable disturbances, the 
acoustic disturbances itself may induce the 
instability waves by interacting with the growing 
boundary layer. Hence to investigate the effects 
of the acoustic and the roughness calculations 
are performed with and without the roughness 
elements. 


Interaction of two-dimensional acoustic 
waves with a boundary layer. 

The model consists of a semi-infinite flat plate 
with a blunt leading edge. The leading edge is 
modeled as a super ellipse of the form 
(x-af y 2 
a 4 b 2 

Here b is the thickness of the plate and in 
accordance with the experiment it is taken as 
.001 inches. The aspect ratio a/b is taken as 10 
hence the blunt leading edge is joined with the 
straight portion of the plate at x=. 01 inch, which 
is at a Reynolds number of 100. 

Figure 15 shows the meanflow density contours 
computed using the WENO code. Figure 14a 
shows the entire domain and figure 14b shows 
the flow field near the leading edge. The leading 
edge shock is located approximately at .0002 
inches, which is 1/5 leading edge thicknesses 
upstream. Computations are also performed with 
an isolated two-dimensional roughness element. 
The roughness is located at x=l in. (. 

( -^Re t = 1000 ) and the height of the roughness 
is h=6/16. Figure 16 (a, b) show the density 
contours including the leading edge shock and 
the Mach wave originating from the roughness. 
The flow field near the roughness is same as that 
computed earlier. The density profiles at x=.01, 

.05, .1, .5, 1., 3. inches. (^e7 = 100, 224,316, 
707, 1000, 1732 ) in the similarity coordinates 
are plotted in figure 17a. In figure 17b, the same 
profiles are plotted in similarity coordinates. The 
Blasius similarity profile is also included for 
comparison. It is seen that very close to the 
leading edge, there exists a strong shock and the 
compression expands over the leading edge and 
becomes weaker away from the nose region. The 
boundary layer profiles slowly approach the 
similarity profiles and even at X=3 in., the outer 
part of the profiles did not merge with the 
similarity profile. 

After the mean flow is computed two- 
dimensional acoustic disturbances are 
introduced at the outer boundaries. The 
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computations are performed at different 
incidence angles 0, 1.5, 10, and of 30 degrees 
and the non-dimensional frequency is F=1.25E- 
4. The amplitude of the acoustic wave is 
p ac lp„ = . 001. 

Figures 18-21 show the results for the incident 
angle of 0 degrees with the roughness element. 
Figure 18a shows the contours of the density 
fluctuations in the whole computational domain 
and figure 18b depicts the results inside the 
boundary layer. Figure 19 shows the maximum 
density fluctuations along the X direction with 
and without the roughness. First observation is 
that the roughness did not introduce any 
difference in the generation of disturbances 
inside the boundary layer compared to that 
generated by the acoustic disturbances alone. 
Hence to the first order, the disturbances inside 
the boundary layer are due to the direct 
interaction of the acoustic disturbances with the 
boundary layer. In figure 20a, the distribution of 
the amplitude of disturbances obtained from the 
simulation is plotted at different axial locations 

x= .5, .1, 3. inches, (-^Re^ = 707, 1000, 1732) 
are shown for zero incident angles. The 
eigenfunctions obtained from linear stability and 
from the continuous spectrum are also included 
for comparison. The question is whether the 
disturbances inside the boundary layers are the 
eigenmodes or the disturbances forced by the 
acoustic waves. The eigenfunctions distribution 
obtained from the linear stability results and 
from the continuous spectrum are similar in this 
region and further analysis is needed before any 
conclusion is made. An important observation is 
that the disturbances are generated very close to 
the leading edge region and after that they 
continues to grow. Figure 21 shows the 
maximum density fluctuations for different 
incident angles 0, 10 and 30 degrees. It is seen 
that the disturbances generated inside the 
boundary layer at 0 and 10 degrees do not 
depend on the direction of the acoustic wave. 
This may be due to the symmetry of the 
problem. The acoustic waves are coming from 
the top and the bottom sides of the plate and 
since the disturbances are generated near the 
leading edge of the plate, the incidence angle 
does not make any noticeable difference in the 


amplitude of the disturbances. The disturbances 
generated by the acoustic waves at 30 degrees 
incident angle initially follows the amplitude for 
the smaller incident angles and decreases 
downstream before it starts to grow again. In 
figure 20b , the eigenfunctions obtained from the 
simulation, linear stability and the continuous 
spectrum are plotted. 


Discussion and Conclusions 

The receptivity process induced by the 
interaction of acoustic disturbances at different 
incident angles over a blunted flat plate with and 
without an isolated roughness is numerically 
investigated. 

The steady field induced by the two and three- 
dimensional roughness is computed. The flow 
field induced by two-dimensional roughness 
exhibits different characteristics depending on 
the roughness heights. A smooth compression 
exists in the front of the roughness, a sharp and 
strong expansion occurs above the peak of the 
roughness and a strong compression takes place 
at the back of the roughness. At small roughness 
heights the flow passes smoothly over the 
roughness, at moderate heights the flow 
separates downstream of the roughness and at 
larger roughness heights the flow separates 
upstream and downstream of the roughness. The 
velocity defects is larger upstream of the 
roughness than that in the downstream. As 
expected the flow field induced by the three- 
dimensional roughness is weaker compared to 
the two-dimensional roughness. 

The interaction of the acoustic waves with a 
blunted flat plate investigation showed that the 
instability waves are generated very close to the 
leading edge. The amplitude of the disturbances 
generated is about five times of the free stream 
acoustic disturbances. It is also revealed that the 
small incident angle of the acoustic wave which 
are impinging the plate symmetrically from the 
above and below the plate does not make much 
difference in the amplitudes compared to that 
produced by a wave which is parallel to the 
plate. Computations have to be continued further 
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downstream and have to be performed with 
asymmetric acoustic disturbances to see the 
effects of incident angles in generating the 
unstable waves. It is also found that isolated 
roughness do not contribute much in generating 
unstable disturbances compared to that is 
generated by acoustic disturbances. Further 
computations have to be done with distributed 
roughness with different length scales to verify 
this conclusion. 
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Figure 6. Variation of the density at constant heights Y/d 0 = 1 and 2 
for different roughness heights h/5 0 =1/16, 1/8, 1/4, 1/2. 
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Figure 7. Modification to the U velocity profiles in the upstream 

and downstream locations due to the roughness elememt 
for roughness height h/8, =1/16. 
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Figure 8. U velocity profiles (a) and the modification to the Blasius 

velocity profiles (b) in the upstream and downstream locations 
due to the roughness element of height h/5 0 =1/2. 
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Figur 9. Contours of the normal velocity (V) and the streamline patterns tor flow over 
Gaussian three-dimensional rouehness elements of heights h/5 = 1/16, 1/8. 
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Figur 10. Contours of the spanwise velocity in the Z-Y plane at dilferent sections 

(X-X 0 )/8=-2, 0, 7.5, 190 for flow over a three-dimensional roughness element 
of height h/8=l/16. 


American Institute of Aeronautics and Astronautics 















AIAA 2003-3589 


Incident wave from above y Reflected wave 


\ 

(Pac<0) 


A, (Pac> 

(Pac>0) 





/ y 



Incident wave from below 




Pac= Pamp exp(ia ac x+ip ac y+ie ac z-icot) 


Figure 11. Schematic for the incident and the reflected 
acoustic field. 



Figure 12. Variation of the wavenumber with the 

Reynolds number for the instability waves 
with P=0 and .025. 




Figure 13. Variation of the reflection coefficient with the inclination angle 
at different Reynolds number for P=.0 and .025. 



Figure 14. Density perturbations obtained by solving 
the continuous spectrum problem. 
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(a) 


(b) 



X (in.) X (in.) 

Figur 18. Contours of the unsteady density fluetuations due to the interaction 
of an acoustic wave with a Hat plate with a blunted leading edge 
and a roughness. F=1.25E-4, incident angle 0.0 deg, M=3.5. 



Figur 19. Maximum density tluetuation with and without 
roughness. F=1 .25E-4, incident angle O.Odeg. 



Figure 20a. Eigenfunctions obtained from simulation and 
linear stability at X=.5, 1., 3 inches. 0 deg. 




Figure 20b. Eigenfunctions obtained from simulation and 
linear stability at X= 1., 3 , 3.9 inches. 30 deg. 


Figure 2 1 . Maximum density fluctuations at different 
incident angles 0. 1 0 and 30 degrees. 
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